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We study cosmologies based on low-energy effective string theory with higher-order string corrections to 
a tree-level action and with a modulus scalar field (dilaton or compactification modulus). In the presence 
of such corrections it is possible to construct nonsingular cosmological solutions in the context of Pre-Big- 
Bang and Ekpyrotic universes. We review the construction of nonsingular bouncing solutions and resulting 
density perturbations in Pre-Big-Bang and Ekpyrotic models. We also discuss the effect of higher-order 
string corrections on dark energy universe and show several interesting possibilities of the avoidance of 
future singularities. 



1 Introduction 

String theory has continuously stimulated its application to cosmology in a number of profound ways 
0] 13 O . It is actually very important to test the viability of string theory by extracting cosmological 
implications from it. In particular string cosmology has an exciting possibility to resolve the big-bang 
singularity which plagues in General Relativity. 

The Pre-Big-Bang (PBB) model lEJIU based on the low energy, tree-level string effective action is one 
of the first attempts to apply string theory to cosmology. In this scenario there exist two disconnected 
branches, one of which corresponds to the dilaton-driven inflationary stage and another of which is the 
Friedmann branch with a decreasing curvature. Then string corrections to the effective action can be im- 
portant around the high-curvature regime where the branch change occurs. Ekpyrotic/Cyclic cosmologies 
[6. ^ have a similarity to the PBB scenario in the sense that the description in terms of the tree-level 
effective action breaks down around the collision of two branes in a five dimensional bulk. 

When the universe evolves toward the strongly coupled, high-curvature regime with growing dilaton, it 
is inevitable to implement higher-order string corrections to the tree-level action. Indeed it was found that 
two branches can be smoothly joined by taking into account the dilatonic higher-order string corrections 
in the context of PBB f^" "101 and Ekpyrotic 1 1 1 1 scenarios. In the system where a (compactification) 
modulus field is dynamically important rather than the dilaton, Antoniadis et al. showed that the big-bang 
singularity can be avoided by including the Gauss-Bonnet (GB) curvature invariant coupled to the modulus 

C3. 

In order to test these string-motivated models from observations, it is important to investigate spectra 
of density perturbations and to compare them with temperature anisotropics in Cosmic Microwave Back- 
ground (CMB). For example, inflationary cosmology generically predicts nearly scale-invariant spectra 
of density perturbations. This prediction agrees well with the recent observations in CMB anisotropics 
1131 . While inflationary cosmology is based upon the potential energy of a scalar field with a slow-roll 
evolution, the kinematic energy of dilaton or modulus field dominates in PBB and Ekpyrotic/Cyclic cos- 
mologies. Hence it is expected that the spectrum of density perturbations in the latter case is different 
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from the prediction in inflationary cosmology. We shall address the problem of density perturbations gen- 
erated in PBB and Ekpyrotic/Cyclic cosmologies by using nonsingular bouncing solutions obtained by 
including second-order string corrections. We note that the effect of string corrections can be important in 
constructing inflation models, see Refs. 1 14 1 for such possibilities. 

The effect of such string corrections can be also important in the context of dark energy. From recent 
observations the equation of state (EOS) parameter w of dark energy lies in a narrow strip around w = —1 
quite likely being below of this value 1151 1161 (see Refs. 1171 for reviews of dark energy). The region 
where the EOS parameter w is less than —1 is referred as a phantom (ghost) dark energy universe [TFl. 
The phantom dominated universe ends up with a finite-time future singularity called Big Rip or Cosmic 
Doomsday 1191 1201 |2T1 . The Big Rip singularity is characterized by divergent behavior of the energy and 
curvature invariants at Big Rip time. Hence it is natural to account for higher-order curvature corrections 
in the presence of dark energy 122. .23. .24. ..25., .26. .27. ..28. ,29, .30, .31. .32] . In fact it is possible to avoid 
or moderate the Big Rip singularity when such corrections are present 1221 l24l 1251 and showed that it is 
possible to avoid or moderate the Big Rip singularity. We shall review cosmological solutions in second- 
order string gravity in the context of dark energy and consider the avoidance of future singularities. 

In what follows the effect of string corrections will be reviewed in two separate sections-(i) PBB/Ekpyrotic 
cosmologies (sec. 2) and (ii) dark energy universe (sec. 3). It is interesting to note that such corrections can 
play important roles for both past and future singularities. 

2 Pre-Big-Bang and Ekpyrotic cosmologies 

The PBB scenario is based upon low-energy, tree-level effective string theory using toroidal compactifica- 
tions j4l|5l. The string effective action in four dimensions is given by 



where is a dilaton field that controls the string coupUng parameter, — e"^, and \/—g is the determinant 
of metric g^^. We neglect here additional modulus fields corresponding to the size and shape of the internal 
space of extra dimensions. The potential Vs{4>) for the dilaton vanishes in the perturbative string effective 
action. The Lagrangian £c corresponds to higher-order string corrections which we will present later, 
whereas £,„ is the Lagrangian of additional matter fields (e.g., fluids, kinetic components, axion etc.). In 
this section we do not consider the contribution of £„i, but in Sec.|3lwe will account for it as a barotropic 
fluid. The above action is so-called the "string frame" action in which the dilaton is coupled to a scalar 
curvature, R. 

The dilaton cj) starts out from a weakly coupled regime (g^ ^ 1) and evolves toward a strongly coupled 
region (gs > 1). The Hubble parameter grows during this stage. This "superinflation" is driven by a 
kinematic energy of the dilaton field, which is so-called a PBB branch. There exists another Friedmann 
branch with a decreasing curvature. It is possible to connect the two branches by accounting for higher- 
order string corrections Cc to the tree-level action i9l ll0ll33ll34l . This is one of the main topics in this 
review. 

In Ekpyrotic J6| Q ™d Cyclic [8J cosmologies the universe contracts before the bounce because of 
the presence of a negative potential characterizing an attraction force between two parallel branes in an 
extra-dimensional bulk. The collision of two parallel branes signals the beginning of the hot, expanding, 
big bang of standard cosmology. After the brane collision the universe connects to a standard Friedmann 
branch as in the case of PBB cosmology. The origin of large-scale structure is supposed to be generated by 
quantum fluctuations of a field (p characterizing the separation of a bulk brane. It is important to construct 
nonsingular bouncing cosmological solutions in order to make concrete prediction of the power spectrum 
generated in Ekpyrotic/CycUc cosmologies. This is actually possible by accounting for higher-order string 
corrections as in the PBB case fl 11 . 
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The PBB model has a similarity to Ekpyrotic/Cyclic cosmologies in a sense that the universe exhibits a 
bounce in "Einstein frame". Making a conformal transformation 

9tJ^u = e^'^Qt^u , (2) 
the action in Einstein frame is given by 

'1 A 1 



Se ^ d xy/-g 



-R--{V4>f-VE{^) 

where Ve(0) = e'^Kg(0). Introducing arescaled field — ±(/)/\/2, the action yields 

'1 ~ 1 



Se — I d^x\/~g 



-R--{vpf-VE{m) + 



(3) 



(4) 



This is the action for an ordinary scalar field ip with potential Ve- Hence it can be used to describe both 
the PBB model in Einstein frame, as well as the ekpyrotic scenario 1 35 1. 

In the original version of the Ekpyrotic scenario ||6l, the Einstein frame is used where the coupling to the 
Ricci curvature is fixed, and the field describes the separation of a bulk brane from our four-dimensional 
orbifold fixed plane. In the case of the second version of the Ekpyrotic scenario [7J and in the cyclic 
scenario (Sj, is the modulus field denoting the size of the orbifold (the separation of the two orbifold 
fixed planes). 

The Ekpyrotic scenario is described by a negative exponential potential (6| 

Fi5 = -yoexp(^-y^<^^ , (5) 

with < p <C 1. The branes are initially widely separated but are approaching each other, which means 
that If begins near +oo and is decreasing toward ip — 0. In the PBB scenario the dilaton starts to evolve 
from a weakly coupled regime with increasing from — oo. If we want the potential (|5} to describe 
a modified PBB scenario with a dilaton potential which is important when but negligible for 

(f) — cx), we have to use the relation ip — — 0/\/2 between the field p in the ekpyrotic case and the 
dilaton </> in the PBB case. 

In the flat Friedmann-Robertson-Walker (FRW) metric ds^ = —dt\ + a\Ax\ in Einstein frame, the 
background equations with £c = are given by 

2HI^]^P^ + Ve{p), (6) 

(p + iHEP^ + -7^ = , (7) 
ap 

where a dot represents a derivative with respect to and He = ciE/aE- Here the subscript "£"' denotes 
the quantities in the Einstein frame. The exponential potential (|3 has the following exact solution L36J 

a.ocM^ He = ^, Ve^-'-^^, P^^. (8) 

tE tE 

The solution for tE < describes the contracting universe prior to the collision of branes. The Ekpyrotic 
scenario corresponds to a slow contraction with < p <C 1. From Eq. (|8ji the potential vanishes for 
p = 1/3, which corresponds to the PBB scenario. 

In string frame the scale factor as and the cosmic time ts are related with those in Einstein frame via 
flie relation dts = e-f/^dts and as = C^/^as- Then we find 111II35I 

asoc{~ts)-^', <l>^-^^\n[-il-^)ts] . (9) 
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This illustrates a super-inflationary solution with growing dilaton. Hence the contraction in Einstein frame 
corresponds to the superinflation driven by a kinematic energy of the field (p. We note that there exists 
another branch of an accelerated contraction {as ex. {—ts)^) t3J, but this is out of our interest. 

The above solution needs to be regularized around ts = (or = 0) in order to connect to the 
Friedmann branch after the bounce. In the context of PBB cosmology, it was realized that higher-order 
string corrections (defined in the string frame) to the action induced by inverse string tension and coupling 
constant corrections can yield a nonsingular background cosmology. A possible set of corrections include 
terms of the form I9lll0l 



£c = -^a'AC(0) [cRIb + di^cj,)^] , (10) 

where a' is a string expansion parameter, C(0) is a general function of (p and Rq^ = R^ — 4_R'"'i?^i, + 
R"^^"^ Raf3fj.u is the Gauss-Bonnet (GB) term. A is an additional parameter which depends on the types of 
string theories: A = —1/4, —1/8 and correspond to bosonic, heterotic and superstrings, respectively. If 
we require that the full action agrees with the three-graviton scattering amplitude, the coefficients c'^s are 
fixed to be c = -1 andd = 1 with = -e""^ |37|. 

The corrections Cc are the sum of the tree-level a' corrections and the quantum n-loop corrections 
(n = 1, 2, 3, • • • ), with the function ^{(f)) given by 

c(0) = -^c„e("-l)^ (11) 

where C„ (ri > 1) are coefficients of n-loop corrections, with Co — 1. There exist regular cosmological 
solutions in the presence of tree-level and one-loop corrections, but this is not realistic in the sense that the 
Hubble rate in Einstein frame continues to increase after the bounce (see Fig. 1 of Ref. 1 11 1). Nonsingular 
bouncing solutions that connect to a Friedmann branch can be obtained by accounting for the corrections 
up to two-loop with a negative coefficient (C2 < 0). 

It was shown in Ref. 1 11 1 that nonsingular bouncing solutions exist in Einstein frame even in the pres- 
ence of a negative exponential potential. When p <C 1 the potential is vanishingly small for ip ^ 1, in 
which case the dynamics of the system is practically the same as that of the zero potential discussed in 
Ref. 1 10|. In this case the dilaton starts out from the low-curvature regime \(j>\ ^ 1, which is followed 
by the string phase with linearly growing dilaton and nearly constant Hubble parameter. During the string 
phase one has \^ 

as a i-Vs)~^, = --^ln(-7/5) + const , (12) 

where (j)f ~ 1.40 and Hf ~ 0.62. In the Einstein frame this corresponds to a contracting Universe with 
aE oc . (13) 

On the other hand, we can consider the scenario where the negative Ekpyrotic potential dominates 
initially but the higher-order correction becomes important when two branes approach sufficiently. By 
including the correction terms of Cc only for < 1, we have numerically confirmed that it is possible to 
obtain regular bouncing solutions, see Fig.^ In this case the background solutions are described by Eq. (|8} 
or (|9|i before the higher-order correction terms begin to work. 

The spectrum of scalar perturbations was studied by a number of authors in the cases of PBB cosmology 
|38 39 40| and Ekpyrotic cosmology |41 42 43, 44| (see also Ref. |45|). A perturbed space-time metric 
has the following form for scalar perturbations in an arbitrary gauge |46 1: 

ds^ = -(1 + 2A)dt^ + 2a(i)B jdxMt + a^it)[{l - 2i/>)% + 2E^,j]dx'dx^ , (14) 
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Fig. 1 Evolution of (p and as with c — —1, d — 1, p — 0.1, Ci = 1.0 and C2 = —1.0 x 10"'^. In this case we 
include the correction term Cc only for (p < 1. We choose initial conditions (f> = —15, H — 1.5 x 10~^. Prior to the 
collision of branes at ip = 0, the universe is slowly contracting, which is followed by the bouncing solution through 
higher-order corrections. 



where a comma denotes the usual flat space coordinate derivative. The curvature perturbation, TZ, in the 
comoving gauge is given by l47l 

TT 

n = ij + ^6(t>. (15) 

The power spectrum of TZ is defined by Vn = fc^|7?.p/(27r^) cx k"'^~^, where fc is a comoving wavenum- 
ber. The spectral index n-ji of curvature perturbations generated before the bounce is given by 1451 14 II (see 
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alsoRefs. I42ll43ll44ll: 



2 



(forO<p< 1/3) 



(16) 



n-R, - 1 = 



l~p 
4- 6p 



(for 1/3 < p < 1) . 



(17) 



We see that a scale-invariant spectrum with n-jz = 1 is obtained either as p ^ oo in an expanding 
universe, corresponding to conventional slow-roll inflation, or for p — 2/3 during collapse 1481 1451 . In 
the case of the PBB cosmology (p = 1/3) one has n-ji = 4, which is a highly blue-tilted spectrum. The 
ekpyrotic scenario corresponds to a slow contraction (0 < p <C 1), in which case we have n-jz ~ 3. 

The spectrum ilbi corresponds to the one generated before the bounce. In order to obtain the final power 
spectrum at sufficient late-times in an expanding branch, we need to connect the contracting branch with 
the Friedmann (expanding) one. As we mentioned, the two branches are joined each other by including 
the corrections given in Eq. il 1\ . This then allows the study of the evolution of cosmological perturbations 
without having to use matching prescriptions. The effects of the higher-order string corrections to the 
action on the evolution of fluctuations in the PBB cosmology was investigated numerically in |49 11). It 
was found that the final spectrum of fluctuations is highly blue-tilted (n-n ~ 4) and the result obtained 
is the same as what follows from the analysis using matching conditions between two Einstein Universes 
1.50.. 5 1 1 joined along a constant scalar field hypersurface. 

It was shown in Ref. fTTi that the spectrum of curvature perturbations long after the bounce is given 
as n-jz ~ 3 for < p ^ 1 by numerically solving perturbation equations in a nonsingular background 
regularized by the correction term dlOl i. In particular comoving curvature perturbations are conserved on 
cosmologically relevant scales much larger than the Hubble radius around the bounce, which means that 
the spectrum (I16> can be used in an expanding background long after the bounce. 

The authors in [7 1 showed that the spectrum of the gravitational potential $, generated before the bounce 
is nearly scale-invariant for < p ^ 1, i.e., n$ — 1 = — 2p/(l — p). A number of authors argued 
[41, 42 , 43 ,'441 that this corresponds to the growing mode in the contracting phase but to the decaying mode 
in the expanding phase. Cartier I 52J recently performed detailed numerical analysis using nonsingular 
perturbation equations and found that in the case of the a' -regularized bounce both $ and TZ exhibit the 
highly blue-tilted spectrum ( I16t long after the bounce. It was numerically shown that the dominant mode of 
the gravitational potential is fully converted into the post-bounce decaying mode. Similar conclusions have 
also been reached in investigations of perturbations in other specific non-singular models 1531 . Arguments 
can given that the comoving curvature perturbation is conserved for adiabatic perturbations on large scales 
under very general conditions |41 , 54J. 

Nevertheless we have to caution that these studies are based on non-singular four-dimensional bounce 
models and in the Ekpyrotic/Cyclic model the bounce is only non-singular in a higher-dimensional com- 
pletion of the model |55|. The ability of the ekpyrotic/cyclic model to produce a scale-invariant spectrum 
of curvature perturbations after the bounce relies on this higher-dimensional physics being fundamentaUy 
different from conventional four-dimensional physics, such that the growing mode of $ in the contracting 
phase does not decay after the bounce L56J . 

3 Dark energy 

In the previous section we have discussed the role of higher-order string corrections to the tree-level action 
in the context of early universe. Recent observations show that the present universe is dominated by 
dark energy responsible for an accelerated expansion. When the universe is dominated by a phantom 
matter (w < —1), this leads to the growth of the energy and curvature invariants. In such a circumstance 
higher-order string corrections may be important when the energy scale grows up to a Planck one. We are 
interested in the effect of such corrections around the Big Rip singularity. In fact the Big Rip singularity 
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can be avoided in the presence of such corrections as we will see in this section. We shall also derive 
cosmological solutions for an effective string Lagrangian together with a barotropic perfect fluid. 
Our starting action is the generalization of Q: 



R) - \uo{m^? - ^(0) + a^w> + ^L^^ 



(18) 



where is a scalar field corresponding either to the dilaton or to another modulus, and / is a generic 
function of the scalar field and the Ricci scalar R. oj, ^ and V are functions of cj). In this section we do not 
consider the cosmological dynamics in the presence of the field potential V. Cm'^ is the Lagrangian of a 
barotropic perfect fluid with energy density p and pressure density p. We assume that the barotropic index, 
w = p/ p, is a constant. In general the fluid can be coupled to the scalar field (j). The a'-order quantum 
corrections are encoded in the term 

C^f = aiRa^p^.R"!'^'' + a2R^,,R'"' + a^^R^ + a^iVclyf , (19) 

where a-i are coefficients depending on the string model one is considering. We are most interested in the 
Gauss-Bonnet parametrization (ai = a-^ — 1, a2 — —4 and — —1) discussed in the previous section, 
but we keep the coefficients general in deriving basic equations. 

For a flat FRW background with a scale factor a, we obtain the Friedmann equation 1491 1251 

d{d - 1)FH^ = RF-f- 2dHF + 2{p + + CPc), (20) 

where d = D — 1 and 

p^ = {lu/2)^^ + V, (21) 

Pc = 3a4^'^ - d[4:CiEH^ + {d-3)ciH'^ + C2{2EHH + 2dH'^H - H"^ + 2HH)]. (22) 

Here F = df/dR, S = ^/^, and 

Ci = 2ai+da2 + d{d+l)a3, C2 = iai + {d + l)a2 + Ma^. (23) 

In four dimensions (d ~ 3), the coefficients read ci = 2ai + 3a2 + 12a3 and C2 — 4(ai + a2 + 803), 
while the i?^ term in pc vanishes. At low energy it was shown that the unique higher-order gravitational 
Lagrangian giving a theory without ghosts is the Gauss-Bonnet one (oi — — 1, a2 = —4, 04 — —1). 
In this case, C2 vanishes identically while ci = 2 + d{d — 3). With fixed dilaton coupling (S = 0) equation 
not reduces to the standard Friedmann equation in four dimensions, in agreement with the fact that the 
GB term is topological when d = 3. In three dimensions (d = 2), the GB higher-derivative contribution 
vanishes identically except for the 0"* term. 

The continuity equation for the dark energy fluid contains a source term given by the coupling be- 
tween this fluid and the scalar field (j). We choose the co variant coupling considered in L57J : 6Sm/54) = 
—y/^ Q{(t>)p, where Sm — J d^Xy/^Cm'' and is an unknown function which we shall set to a 
constant later. In synchronous gauge we have 

p + dHp{l + w) ^ Q^p, (24) 

while the equation of motion for the field cj) is 

u{4> + dH4>) + V' +'^aiC<P'^{34> + dH^ + E4>)+ (^uj^-^^ = -Qp, (25) 

where the Lagrangian of the quantum correction is written as 

d't'^ ^d\{d+ l)ciH* + AciH^H + C2H^] + ai4>'^. (26) 
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Equations ( l20> - (l25t are the master equations of the physical system under study. We note that the ij)'^ term 
can be important even in the absence of curvature corrections in a dilatonic ghost condensate model 1591 . 
The massless dilaton discussed in the previous section corresponds to 

F = -oj = e-^, V^Q, e=^e-^. (27) 

The full contribution of rt-loop corrections is given by Eq. (^3- In this work we shall take only the tree- 
level term ( I27> into account. 

Generally moduli fields appear whenever a submanifold of the target spacetime is compactified with 
compactification radii described by the expectation values of the moduli themselves. In the case of a single 
modulus (one common characteristic length) and heterotic string (A = 1/8), the four-dimensional action 
corresponds to L58il 



F = l, uj = 3/2, a4 = 0, ^ = ln[2e*r/''(ie'^)], (28) 

where 77 is the Dedekind function and (5 is a constant proportional to the 4D trace anomaly. 5 depends on 
the number of chiral, vector, and spin-3/2 massless supermultiplets of the N = 2 sector of the theory. In 
general it can be either positive or negative, but it is positive for the theories in which not too many vector 
bosons are present. Again the scalar field corresponds to a flat direction in the space of nonequivalent 
vacua and = 0. At large the last equation can be approximated as 

C« Co cosh,/), Co=^, (29) 

which we shall use instead of the exact expression. In fact it was shown in Ref. |,60 1 that this approximation 
gives results very close to those of the exact case. 

3.1 Modulus driven solution 

In Ref. [25] cosmological solutions based on the action ( I18> without a potential {V = 0) were discussed 
in details for three cases-(i) fixed scalar field (</) = 0), (ii) linear dilaton (0 = const), and (iii) logarithmic 
modulus (0 cx 1/t). In the case (i) we obtain geometrical inflationary solutions only for £) ^ 4. In the case 
(ii) pure de-Sitter solutions exist in string frame, but this corresponds to a contracting universe in Einstein 
frame. These solutions are not realistic when we apply to dark energy. In what follows we shall focus on 
cosmological solutions in the case (iii) with a fixed dilaton. 
Introducing the following new variables 

X = H , y = H , u = (j), V = 4> , z = p , (30) 

the equations of motion for the modulus action corresponding to Eqs. J20t -( l25t read 



y 



i^y, (31) 

2z + 3v'^/2 - d{d - l)a;2 - 8dci£,x^ - 2d{d - 8)01^2;'* - 2dc2^y(2Ex + 2dx'^ - y) 

Adc2^x ' 

(32) 

'^Ciid + l)cix^ + ^cix^V + C2y2] _ dxv - ^z, (33) 



V 

z = [-dx{l + w) + Qv] z. (34) 



While only derivatives of ^ appear in the equations of motion for d — "i and C2 — Q (GB case), there are 
non-vanishing contributions of ^ itself for general coefficients c^. When C2 — 0, the equations of motion 
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for X and v read 



z + 3vv/2 - dcix^ [4g +{d~ 3)ix] 
dx[id - 1) + 12ciix + 4(d - 3)ciS,x^] ' 



V 



cie'^[(d+l)x2+4i]-(i-- 



X 2(5 z 



3 v2' 



3 w 

while the Friedmann equation is 

d(d - 1) 4 - ? - 24 + 2dci^4 [42 + (d - 3)a;] = 0. 



(35) 
(36) 

(37) 



In addition ci can be set equal to 1 and absorbed in the definition of ^0, so that the coefficient C2 is the only 
free parameter of the higher-order Lagrangian. 

We search for future asymptotic solution of the form 



UJl 



zot^'^^ exp 



d{l + w)uji 
/3+1 



z ~ zot", 13 =-1, 
where the barotropic index w is constant and 

a = Qu-i — d{\ + vS)lo\. 



(38) 

(39) 
(40) 

(41) 



We define 6 = (l/2)ci^oe"''j so that in any claim involving the sign of 8 (S) a positive c\ coefficient is 
understood. In order to find a solution in the limit t +00, one has to match the exponents of t to get 
algebraic equations in the parameters /3, Wi, C2 and (5. 

In Ref. L58..25J the following four solutions were found. We can obtain a number of analytic solutions 
depending on the regimes we are in: 

1. A low-curvature regime in which ^ terms are subdominant at late times. 

2. An intermediate regime where some terms in the equations of motion, either coupled to ^ or not, are 
damped. 

3. A high-curvature regime in which ^ terms dominate. 

4. A solution of the form j38>-j40> for the full equations of motion. 

Below we summarize the properties of each solution. 

1 . A low-curvature regime 

In this regime the solution is given by 

= c^2<2, (42) 

with the constraints 



1 IQz^t 



Q + 2 



3dw2 



3cj2 = 2d{d - - 4zoi 



Q + 2 



a < -2. 



(43) 
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2. An intermediate regime 

In this regime the solution exists only for C2 — and is given by 

U2^5, Q<-2/5, c^? = (15 - 2Qzoi^^+2) , (44) 

Ibao 

for a non-vanishing fluid. The condition, ^ > 0, is required in order to obtain an expanding solution 
characterized by a{t) ~ gq exp(— wi/i). This solution reaches Minkowski spacetime asymptotically. 
If the fluid decays, then one recovers the Coo solution of Ref. 115 81 with d = 3 and z = Q. 

3. A high-curvature regime 

In this regime the solution is given by 

/3=-l, W2 > 2, a<w2-4, (45) 

together with constraint equations 

6dul[{d + 1)luI - Auji + C2] - Qzot"-'^''+'^ = 0, (46) 
d5ujl[{3 - d)ujl - 4w2Wi + C2(2du;i + 2^2 - 3)] + zoi""""+^ = 0. (47) 

In the GB case (c2 = 0) the solution corresponding to a decaying fluid (a < 0^2 — 4) is 

4 3-d 

"^ = dTT' "^ = dTT' ^^^^ 

which contradicts the condition ( 145 1 in any dimension. Hence in the GB scenario with 0^2 7^ 4 only 
the marginal case a = (jJ2 ^ 4 is allowed. 

4. An exact solution 

An exact solution which is valid at all times is 

/? = -!, L^2 = 2, (49) 

together with the constraints on ui: 

2ddujf [{d + 1)^1 - 4wi + C2] - 6duji + 6 - QzQt"+'^ = 0, (50) 
2ddLol[id - 3)luI + 8wi - C2(2duji + 1)] + d{d - 1)luI - 6 - 2zot°'+'^ = 0, (51) 

and Eq. ( l43t . 

The low-curvature solution J42t and the Minkowski solution J44t can be joined each other if the coupling 
constant S given in Eq. ( I28> is negative 158, 61J. The exact solution ( I49> is found to be unstable in numerical 
simulations of Ref. |25|. In the asymptotic future the solutions tend to approach the low-curvature one 
given by Eq. i42\ rather than the others, irrespective of the sign of the modulus-to-curvature coupling S. 
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Solutions t ^ oo 




^1 


Stability 


Low curvature 




Yes 


High curvature, C2 — 




A > 


High curvature, C2 7^ 0, A = 






High curvature, C2 ^ 0, A 7^ 


tO^, C2 < -8 


w+, A > 


Exact, C2 = 






Exact, C2 7^ 


Yes 





Table 1 Constraints on modulus solutions in the asymptotic future for the Hubble parameter H = uit ^, vanishing 
fluid, and d — 3. Blank entries are excluded by experiments or numerical analysis. 



3.2 Constraints from the recent universe 

We compare the observational constraints on uji for the recent evolution of the universe with the modulus 
solutions found in the previous section (d — 3). The situation we study is the case in which a perfect 
fluid is vanishing asymptotically. The results are summarized in Table^at the 68% confidence level. We 
also address the cases in the presence of a cosmological constant A. Note that the solution i44i in an 
intermediate regime is discarded. 

The logarithmic modulus solution with the GB parametrization and no extra fluid does not provide a 
viable cosmological evolution in the current universe. In the next subsection, however, we will see that the 
GB case in the presence of dark energy fluid may exhibit interesting features for the future evolution of the 
universe. The low redshift constraint on C2 for the high curvature solution (A 7^ 0) can be relaxed up to 
C2 < — 1 at the 99% confidence level. Hence we have shown that there are models which can in principle 
explain the current acceleration without using the dark energy fluid. 

The situation becomes more complicated in the presence of a barotropic fluid. The low-curvature so- 
lution can describe the very recent universe if Q is negative and non- vanishing. The other cases crucially 
depend upon the interplay between all the theoretical parameters. 

3.3 Dark energy universe with modulus gravity 

In the universe dominated by a phantom fluid (w < —1), the energy density of the universe continues to 
grow and the Hubble rate eventually exhibits a divergence at finite time (Big Rip). Then the effect of higher- 
order curvature corrections can become important when the energy density grows up to the Planck scale. 
In fact it was shown that quantum curvature corrections coming from conformal anomaly can moderate the 
future singularities 1 63 22 1. 

We would like to consider the effect of a' quantum corrections when the curvature of the universe 
increases in the presence of a phantom fluid. We shall concentrate on the modulus case with given by 
Eq. i29\ . Our main interest is the cosmological evolution in four dimensions (d — 3) in the presence of a 
GB term. The dilaton is assumed to be fixed, so that there are no long-range forces to take into account 
except gravity. 

From the discussion in subsection 13. II the growth of the barotropic fluid is weaker than that of the 
Hubble rate when the condition 

a = Quj2 - d{l + w)uji < -2 , (52) 

is satisfied. This condition is not achieved for a phantom fluid when the coupling Q between the fluid and 
the field (j) is absent {Q = 0). 

In Ref. P251 the equations of motion were solved numerically by varying initial conditions of H, (j) and 
p. When (5 < 0, we numerically found that the solutions approach a Big Rip singularity for Q — and w < 
— 1 (see Fig. The condition ( I43t can be satisfied for negative Q provided that uj2 is positive. In Fig.|2]we 
plot the evolution of H and p for Q = —5 and w = —1.1. In this case p decreases faster than which 
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0.10 
0.010 
0.0010 
0.00010 
1 0" ' 
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t 

Fig. 2 Evolution of H and p with = —2, w — —1.1 for (sl) Q — and (b) Q = —5. We choose initial conditions 
as Hi = 0.2, (f>i — 2.0 and pi — 0.1. The curves (al) and (bl) represent the evolution of if for Q = and Q = —5, 
respectively, while the curves (a2) and (b2) show the evolution of p for corresponding Q. 



means that the energy density of the fluid eventually becomes negligible relative to that of the modulus. 
Hence the universe approaches the low-curvature solution given by Eq. (142 > at late times, thereby showing 
the avoidance of Big Rip singularity even for w < —1. By substituting the asymptotic values cui = 1/3 
and UJ2 = 2/3 in equation ( I43> . the condition for decaying fluid reads Q < 3(w — l)/2 = —3.15. We 
checked that the Big Rip singularity can be avoided in a wide range of the parameter space for negative Q. 
These results do not change even for smaller values of 5 such as \S\ = 0(1) (corresponding to |^o| = 0.1). 

When Q is positive, the condition ( I43> is not satisfied for > 0. However numerical calculations 
show that (j) becomes negative even if > initially. We found that the system approaches the low-energy 
regime characterized by wi = 1/3 and llI2 = — 2/3. Since u}2 < 0, the Big Rip singularity can be avoided 
even for positive Q. In fact we numerically checked that the Hubble rate continues to decrease provided 
that the condition ( I43t is satisfied in the asymptotic regime. 

When 5 > 0, there is another interesting situation in which the Hubble rate decreases in spite of the 
increase of the energy density of the fluid. This corresponds to the solution in the high-curvature regime in 
which the growing energy density p can balance with the GB term (p « 24iJ^^ in the Friedmann equation). 
In this case the Big Rip does not appear even when w < —I and Q — 0. Thus the GB corrections provide 
us several interesting possibilities to avoid the Big Rip singularity. 

3.4 Other corrections 

We should mention several attempts to apply curvature corrections to dark energy. Nojiri et al. 1231 
studied the case of the tree-level dilatonic-type coupling (I27> in the presence of an exponential potential 
and showed that the scalar-Gauss-Bonnet coupling acts against the occurrence of the Big Rip singularity. 
In Ref. 1241 the evolution of phantom dark energy universe was studied with fixed dilaton and modulus 
by taking into account string corrections up to the 4-th order It was found that the universe reaches a Big 
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Crunch singularity instead of a Big Rip one for type II strings, while the Big Rip singularity is not avoided 
for heterotic and bosonic strings. 

Apart from string corrections, a number of authors 1621 1631 1221 l64l 1651 studied the effect of quantum 
backreactions of conformal matter around several singularities which can appear in future. They usually 
contain second-order curvature corrections such as the Gauss-Bonnet term and the square of a Weyl tensor. 
In Ref. L22J it was shown that quantum corrections coming from conformal anomaly can be important 
when the curvature of the universe grows large, which typically moderates future singularities. Finally we 
note that loop quantum cosmology leads to a modified Friedmann equation when the energy scale grows to 
a Planck one, which typically gives us a regular cosmological evolution without future singularities 1.66J . 

4 Conclusions 

In this article we have discussed cosmological implications of higher-order string corrections to the tree- 
level effective action. In the context of Pre-Big-Bang (PBB) and Ekpyrotic cosmologies regular bouncing 
solutions can be constructed by including such corrections. This allows us to evaluate the spectrum of 
density perturbations long after the bounce. For the correction terms given by Eq. (II 1> we found that 
the spectra of scalar perturbations are highly blue-tilted: n-jz = 4 in the PBB case and n-ji = 3 in the 
Ekpyrotic case. This is different from the nearly scale-invariant spectrum (n-ji ~ 1) observed in CMB 
anisotropics. As long as nonsingular bouncing solutions are constructed by using the correction terms 
presented in this paper, we need another scalar field (e.g., curvaton |23) to generate nearly scale-invariant 
density perturbations. 

We also applied second-order string corrections to dark energy. In particular we reviewed several cos- 
mological solutions in the presence of modulus-type corrections with a fixed dilaton. In the asymptotic 
future the solutions tend to approach the low-curvature one given by equation (I42t rather than the others, 
irrespective of the sign of the modulus-to-curvature coupling S. We placed constraints on the viability of 
modulus-driven solutions using the current observational data. The Gauss-Bonnet parametrization is ex- 
cluded in any of the above mentioned regimes when a barotropic fluid is vanishing; see tabled In the pres- 
ence of a phantom dark energy fluid we discussed the effect of the modulus coupling with Gauss-Bonnet 
curvature invariants. It is possible to consider a situation in which the energy density of the fluid decays 
when the coupling Q between the field and the phantom fluid (w < —1) is present. We showed that the 
Big Rip singularity can be avoided for the coupling Q which satisfies the condition Qijj2 — <i{l+w)u!i < — 2 
asymptotically. This is actually achieved irrespective of the sign of Q and the asymptotic solutions are de- 
scribed by the low-curvature one given by Eq. (I42> . We also briefly mentioned the effect of other forms of 
higher-order string corrections to future singularities. 

Thus we showed that string corrections can be important in a number of cosmological situations. We 
hope that the development of string theory will further provide us rich and fruitful implications to cosmol- 
ogy- 
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